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Abstract 

Stability of standing flux-antiflux front in anisotropic layered superconductors is considered. We 
describe two assisting mechanisms destabilizing the standing vortex-antivortex front. There are 
anisotropy of the layered superconductors and the heat, released by the vortex dynamics. We 
present the conditions of the front stability for various anisotropy and heating parameters. We 
predict that even small vortex-antivortex heating can result in front instability. The characteristic 
size of the unstable pattern is estimated. 
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I. INTRODUCTION 



Studies of patterns in the magnetic flux distribution in the type-II superconductors are 
attracting the attention of many research groups [l-4] whose magneto-optical experiments 
demonstrate that nonuniform flux penetration occurs. Patterns with branch-like structures 
lave been found in most of high T c materials, like y_Ba 2 CM307_a;j^] and Bi2Sr 2 CaCu2 Og +x 
6]. The nucleation of dendrite like patterns in MgB 2 films is another example [z, 8|. These 
complex structures consist of alternating low and high vortex density regions and are found in 
a certain temperature window. Likewise, flux penetration in the form of droplets separating 
areas of different densities of vortices has been observed in NbSe2 jsj. Usually the occurrence 
of flux patterns in interfacial growth phenomena can be attributed to a diffusion driven, long- 
wavelength instability of a straight front, similar to the Mullens-Sekerka instability 10| found 



in crystal growth. The nucleation of nonuniform patterns associated with the propagation 
of a flux front into a flux-free sample has been attributed to such an interfacial instability. 
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13] where a higher temperature leads 



This results from a thermomagnetic coupling [7, 8, 
to a higher vortex mobility, enhanced flux flow, and hence a larger heat generation. 

On the other hand, the situation when the vortices interact in a superconductor with 
the flux of the opposite sign is less theoretically studied. This flux configurations arises for 
example when a DC bias current creates vortices and antivortices on the opposite side of 
the superconducting strip [14]. Another example which is now under intensive investigation 
arises upon exposing the previously magnetized sample to the magnetic field of an opposite 
direction. According to the experimental data (2 4, 15, [l6], the boundary between vortices 
and antivortices exhibits many features suggestive of a long wavelength instability. The cause 
of the instability at the boundary between fluxes of opposite sign is still being debated. In 
particular, Fisher at el. 17J, |l8) proposed a non thermomagnetic mechanism of instability 
caused by an in-plane anisotropy of the vortex mobility. This mechanism of instability was 



carefu 
el. 171 



ly reinvestigate by the van Saarloos et el 19] . They confirmed the finding of Fisher et 



18] that standing vortex-antivortex fronts have an instability to a modulated state, 



while the moving fronts were found to be stable for all anisotropies. In fact however, the 
fl— ux ^taMity was experimentally detected in a system with small Q and moderate 
anisotropy [3]. Several years latter this model was improved by an additional assumption 
of a step shape and anisotropy of the voltage current characteristics [18] and explained the 
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experimental result in moderate anisotropic superconductor Y BCI2CU3O7-5 . Unfortunately 
this assumption cannot explain the instability in pure isotropic systems like Nb and MgB2 . 
On the other hand the thermomagnetic mechanism can also be responsible for flux-antiflux 
instability when the flux-antiflux front is heated both by the vortex (antivortex) dynamics 
and by the vortex-antivortex annihilation. 

In this paper we report on the thermomagnetic theory of the flux-antiflux instability in 
anisotropic layered superconductor. 



II. MODEL AND BASIC EQUATIONS 

n 

We start with a model of two-component vortex gas [20(] spatially homogeneous along 
the z axes, which is valid for the experimentally interesting situation of the low magnetic 
field when typical spacing between vortices a® essentially exceeds vortex-vortex (antivortex) 
interaction radius £, and the vortex velocity depends only on the edge screening current 
that is assumed to be homogeneously distributed across the sample. The vortex- vortex 
repulsion, in this case, keeping the number of vortices, cannot play a significant role. One 
must take into account both vortex-antivortex annihilation and heat release accompanying 
this process. We should also take into consideration heat absorption by the sample lattice 
in order to prevent the rise of unlimited temperature. 

The vortex-antivortex annihilation obeys the well-known master equations of the recom- 



bination theory 2l|, [22 ] 



^+V(n + v+)+#n + n_ = 0, (1) 

— + V(n_v_)+£n + n_ = (2) 

g = f Vj v = mod (v+ - v_) (3) 

where n + and n_ are the vortex and antivortex densities, respectively, g is the ratio of 
recombination for vortices and antivortices, £ is the cross section of the annihilation, which 
is of the order of the coherence length of the superconductor, and v± are the opposite 
directed vortex-antivortex velocities, which in the creep regime are strongly temperature 
dependent: 

mod(v ± ) =v± = v± FF exp (- — 



3 



J ~V x B, B = ip Q (n+ - n_) , v ±FF = ll^. (4) 

47T 77c 

Here, [/ is a temperature-dependent pinning potential, <£> is the unit flux, J is the electric 
current, and r\ is the viscosity of the vortices, which in an anisotropic system is a tensor 
with different in plane and across the plane tensor components r) a p. 

(It should be noted that our master (recombination) equations account numbers of the 
topological charges (vortex/antivortex cores). The number of the vortices and antivortices is 
changed when they meet each other at the distance of the order of the coherence length £ , 
rather than the penetration length of the magnetic field A. For example, the same equations 
describe vortex-antivortex annihilation both in the superfluid He4 and in superconducting 
films where the magnetic field is uniform). 

We wish to investigate an anisotropic system where the vortices velocity is not necessarily 
parallel or perpendicular to the layers of the superconductor. Therefor, we shall take the 
vortices velocity in a general form 

Va = ~7«fl ( I" J • ( 5 ) 



a,/3 = x,y 

Here 7 = (p /47ri] when 



. , cos 2 1? + a sin 2 -d cos $ sin $ (1 — a) . 

7 = V = Vo ( I ( 6 ) 

cost? sin $ (1 — a) a cos 2 •& + sin 2 •& 



is the inverse tensor of the vortex viscosity [19], a is the anisotropy parameter of the system 
(0 < a < 1), d is the angle between the x axis and the a — b plain of the layered structure 
(Fig. 1). 



12] one can corn- 



Assuming that the heat diffusion length exceeds the width of the slab 
plete the set of equations Eqs. flTJ and fl2]) by the temperature transfer equation in the 
form 

9T Kt 1 6Q (T-Tp) 

dt C p C p 8t t R ' 1 ' 



where 



d -9- = Wj + W A (8) 
Wj = ^- (n+ + n_) ; W A = Zv^Qo (9) 
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is determined by the energy released both by vortex-antivortex dynamics Wj and by vortex- 
antivortex annihilation Wa- Here Kt is the heat conductivity, C p is the heat capacity, To 
is the coolant temperature, Qq is heat released by annihilation of a single vortex-antivortex 
pair per unit vortex length, and £jj is the characteristic time of temperature relaxation. 

The set of Eqs. (JTJ)-([T]) completed by the boundary conditions describes all features of 
the model. 



III. SPATIAL DISTRIBUTION OF FLUX-ANTIFLUX DENSITIES 

We consider the case of a restricted sample of length D. In this case, the flux-antiflux 
interface in the stationary state is formed due to a balance between flux-antiflux entering 
the sample from the opposite sides and their annihilation in the middle point. 

Introducing new variables 



AL, x/AL -> x, 

2 



n+/n m = N+, n_/n ni 

AL=(n m £T\ t/t Q =t\ t 

b = N+-N_, N = N+ + N. 



2 ' 



(10) 

(11) 



where n m is flux density at the interface point, and AL is the characteristic width of the 
region in which the spatial distributions of vortex and antivortex flux densities overlap, 
forming the interlayer where the vortices of the opposite signs coexist. AL may be estimated 
as AL ~ v/gn m ~ (n m 1 (see Ref. [23]), which is a microscopically large area where the 
total magnetic induction is suppressed. 
One obtain from Eqs. (0Q)-(j2]) 



db 


d 


dt* 


dx 


-{ 




dx \ 


^ dx ) 



,db\ 

dx J 



dN 
dF 

where b is the dimensionless magnetic induction. 
In the stationary state we get for N and b 

dx V dx 



\N 2 ~b 2 \ 



db 



dx 



(12) 
(13) 



(14) 
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Performing the integration in Eq. (j!4|) one obtain 

N<* = _,. (16) 

here J is a constant. 

Substituting A" from Eq. (fT6l) into Eq. f|T5|) and performing the integration we immedi- 
ately obtain the differential equation for b in the form 

I + W 



W + -In 



I-W 



b 3 

1- (17) 



where 



W = 6o §. ( 18 ) 



Flux-ant iflux Interface 



This equation may be solved analytically close to the interface line where bo goes to zero. 
Assuming that the vortices and antivortices appear at the edges of the samples separated 
by the distance D (in dimensionless units) and assuming the following boundary conditions 

u dh » T at n + D 
b — = -I, A/ _ = at x = - — , 

ox 2 

b ^ = I, N 0+ = at x = ^, (19) 

we obtain an asymptotically exact result for magnetic induction and vortices density at the 
interface: Looking for the solution in the vicinity of the ffux-antiflux front (X — > 0) in the 
form 

bo ~ a±x + a2X 3 ; Ao ~ 03 + a 4 x 2 (20) 

one obtains from Eqs. f fT6|) -f fT8|) to the main order 

7-4/3 7 

a x = —J 2 / 3 , a 2 = , a 3 = J 1/3 , a 4 = — (21) 

'180 60 V ; 

It should be noted numerical simulation show that these formulas are valid in a much wider 
region at front (b = 0) and can be considered as an interpolation ones. 



Assuming that the slope of the magnetic induction at the front / 1 is small one obtains 
for the characteristic size of the vortex-antivortex area x c , where vortices and antivortices 
coexist N (x c ) = b (x c ) (see Eq. lfHj) ) (see Fig. 2). 



x r 



j-i/3 x (22) 



The dimensionless vortex velocity at the interface uq = I 2 / 3 is a constant. Returning to 
the dimension variables, we obtain for interface flux velocity 



w n^gexp(-y/T) /2/3 

4"7T 7] 



IV. OVERHEATING INSTABILITY 

Let us consider the stability of the vortex-antivortex interface with respect to small 
deviations from its initial plane shape. For this we shall take the Eqs. ( fl2l) . ( 1T31) . and (JTj) in 
the more general form 

| - V (Ny) = (24) 

V(6v) + (N 2 -b 2 ) \v\ = (25) 

— - kV 2 6 - w A - wj + r (0 - 1) = (26) 

where 

= {N 2 - b 2 ) S A \v\ , wj = NSjv 2 (27) 

are the dimensionless annihilation and Joule heat terms. 

Here Sa = (Qo n m/^T C p ) ; Sj = (pln 2 m /16ii 2 T Cp are the heating parameters, Qo ~ 
<y9g/A 2 where A is the London penetration length. The ratio Sj/Sa ~ A 2 n m . 

It seems at first glance that the direct Joule term caused by vortex (antivortex) motion 
always prevails. Really, for a sharp shape magnetic induction front the vortex-antivortex 
annihilation term (overlapping) which is proportional to the vortex (antivortex) density 
production wa ~ vN + N_ is small while vortex velocity v ~ V6 is large. Therefore the 
direct, Joule term which is proportional both to the sum of the vortex and antivortex 
densities and to the square of the velocity wj = NSjv 2 significantly exceeds the annihilation 
term. However, in our case, when the slope of the magnetic induction profile is small the 
annihilation term becomes essentially important. In this case the overlapping (annihilation) 



term is larger due to deep mutual penetration of vortices and antivortices over the interface 
area (see Fig 2). The vortex velocity in this case is small and it decreases the Joule term 
which is of the order of v 2 . Substituting functions N + ,N_ from the Eqs. (|2T)|) .( l2T|) . (lTTi) into 
Eq.f[2"Tl) one obtains for the wj ~ N (V6 ) 2 — I 5 ^ 3 and wa ~ N^Vbo ~ J 4//3 allowing to 
neglect in our consideration the Joule term which is relatively small wj/wa — I 1 ^ 3 << 1. 
Here the dimensionless velocity v has the form 

(28) 

while k — > t KT/cp(AL) z (here Kd = kt/C p is the diffusion constant), r — > t /tn are the 
dimensionless effective diffusion and relaxation coefficients correspondingly. 



v = exp (-U/T) 



db 

dx 



A. Small fluctuations 



Looking for a solution of the form 

b(x,y,t) =b (x) + i/j (x, y, t) 

N(x,y,t) = N + C(x,y,t), 

Q(x,y,t) = l + 9(x,y,t), 

where ip, ( and 9 are the small perturbations of the form 

/ 



(29) 

(30) 
(31) 



i/j (x,y,t) 
C(x,y,t) 
y 6(x,y,t) J 



» ' $ (x) » 



exp (Xt + iky) , 



(32) 



C(x) 

\e(x) J 

Here ip Ql Co an d are constant amplitudes, A is the rate grow of the perturbation and k 
is the wave vector in the y direction. 

Taking into account the fluctuations of the vortex velocity 

db 



(33) 



where 9/T is the change of velocity due to thermal fluctuations [23[ and 7 = <p Q /4:iir] (see 
Eq. [6])) one obtains 



db dtp db dip 



' xy 



dx lxx dx lxx dx lxy dy 
_ dbo _ dj)_ _ ^0n_ 

Vy + 6v y~ 7 ^ Q x V Qx lyx dx lyy dy 



dx 
dijj 
dx 



dx 
dbo 
dx 



(34) 
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Here 5v x , y are the deviations from the steady state vortex velocity. 

While the vortex velocity in the x direction is higher then in the y direction, we assume 
that \v\ ~ v x . 

We neglect the influence of temperature fluctuations on heat capacity C p and relaxation 
coefficient r because their calculations do not result in essential effects. We also neglect in 
the main order the change of the average temperature in the flux front area. 

Substituting the perturbations in the form [32] into initial set of Eqs. fl24l) - fl26l) . and use 
the stationary solution in the form Eqs. ffT4j) - (TT5"]) (see Appendix I) one obtains from the 
Eq. (j581) for the rate grow 

a 3 + a 2 (r\ + n^ 2 ) + a (r 2 + u 2 k 2 + u 3 k 4 ) + (u 4 k 2 + n 5 fc 4 ) = o, (35) 

where 7^ = 7^ and 

T, = 2 lxx I + r- lxx liS A (36) 

n i = -^yy 1 ^ + K 

T 2 = 2 7xx ./r 

n 2 = 2 (-fly - i xx i yy ) /3 + 2 1xx Ik - l yy I^r + liS A (7 M 7 ni - 7^) 
II '. = ' .i.j 1 ' 1 ' 

04 = 2 (7^ - J XX Jyy) Til + J XX jlyS A i% 

KB = 2 {jl y - -1 XX lyy) K/5 

The roots of these equations are presented in Appendix where the solution A12 are rele- 
vant. 

(We consider only the nonuniform instability, hence solutions with Re A >0 at k = are 
omitted). The onset of the nonuniform along the front instability is determined either by 
the conditions (Re Ai = 0) 

(Re£i) 2 = 4 (Re Ci) (Re A x ) (37) 
ReCi < 0, (Re Ai) < (38) 

giving the contact point at Re Ai — k 2 plane (see figs 3-5) 

* 2 - -S (39) 

or at C\ = (Rei?i > 0,ReAi < 0) resulting in the Mullens- Sekerka instability. 



V. RESULTS 



A. Strong Heating. 



We start with a model case when the heating coefficient S — > oo. In this case the param- 
eters of the dispersion equation [36] have the form 



Ti = -~f xx r*s A ] iii = -i w i* + «; r 2 = 27^/7-; n 2 = i*s A (7^7^ - 7^) ; ( 40 ) 



n 3 = -7 ro ^»«; n 4 = 7xx7xj/^a/*; n 5 = 2 (7^ - 7 SKe 7 w ) k/3; -y/rf - 4r 2 = 7 :e:e / 4/3 5a 



while the equation for Re A 2 function reads (see Appendix II) 

ReA2 = 2i^_%& (41) 
4r 64r 3 v ' 

The rapidly growing mode in this case d 2 Re X^jdk 2 = (maximum velocity of the mostly 
unstable mode) has the wave vector k = 2\/2r / xy S A I and the period of the pattern 
along the front (see fig. 3) 

, f3ir-f S A I , . 

dy = \/2 r ( } 

In the case of moderate and small heating parameter S A , the results are strongly depends 
both on anisotropy and on the relation between other parameters and can be done numeri- 
cally. 



B. Moderate and weak heating. 

In this case the equation for real part of the increment of the instability is determined by 
the equation (see Appendix II) 

Re Ai = Re Ci + k 2 Re B l + k A Re A x ; (43) 

It has been solved for different anisotropy, in-plane diffusion coefficients k and relaxation 
coefficients r.The results are presented in figs. 4-6. In all of the curves at these figures the 
Re Ci < 0, Re B x > and Re A x < 0. 
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The instabilities in all of these cases has the form of contact one rather than the Mullens- 
Sekerka type. The nonuniform structure along the front appears with the period d y = 27r/k c 
where k c is the contact point of the Re Ai with k axis. 

In fig. 4 the Re Ai as a function of k 2 is shown for anisotropic superconductors with 
various in-plane diffusion constant k. There is a critical heating parameter S and critical 
diffusion constant when the instability arises (curve 1), while the system becomes stable 
as the diffusion constants increase (curve 2,3). The relaxation constant r (the coefficient 
of the ballistic heat conductivity) also strongly affect the instability condition. In fig. 5 the 
increment Re A versus k 2 for anisotropic superconductor with different relaxation constant 
r demonstrate that the instability appearing at relatively small constant r (see curve 1) 
disappears as the relaxation parameter grow (curves 2,3). 

Fig. 6 demonstrates that the anisotropy essentially affects the onset of the instability. 
The increment of instability Re A versus k 2 for different anisotropy. Curve 1 for isotropic 
superconductor (a = 0.9, i9 = 7r/4, j xx = ^ yy = 0.545,7^ = 0.055) shows the instability at 
heating coefficient Sa = 0.89 while curve 2 demonstrates the lack of instability at heating 
coefficient Sa = 0.8. The curve 3 exhibits instability for anisotropic superconductor (a = 
0.1, $ = 7r/4, j xx = = 0.55, -y xy = 0.45) and even more small heating coefficient Sa = 0.1. 



VI. CONCLUSIONS 



The standing flux-antiflux interface demonstrates instability due to the heat released by 
the flux-antiflux annihilation. In fact this is a well known Kelvin-Helmholtz (KH) instab ility 



appearing when different layers of liquid move with the opposite directed velocities 25]. 
In our case, however, vortex and antivortex "liquids" are moving as it is shown in Fig. 7. 
The heat released by the annihilation enhances the vortex/antivortex velocities resulting in 
turbulence instability at the flux-antiflux interface. The rate grow dependence on the wave 
vector directed along the front showing the instability is presented in Figs. 3-6 for different 
heating parameters and anisotropy of the superconducting materials. The characteristic size 
of the unstable pattern is determined either by the rapidly growing mode of the real part of 
the rate grow d^jdk (for strong heating parameter) (see fig. 3) or at the contact mode for 
moderate and small heating (figs. 4-6). 

The theory predicts stability of the flux-antiflux front for any physical parameter of the 
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system (see Eqs. A12 from Appendix II). The physical reason for the instability is the result 
of growing temperature gradients along the front when vortices are moving with different 
velocities. The velocity of the flux flow vortices is very high and more rapid parts of the 
front can break it during the time of the instability. In particular for materials with typical 
parameters 

B = 2000G, r] = 5 • 10~ 5 CGSE, v F = 10 7 cm/ sec, I = l(T 8 cm, f = 10~ 6 cm (44) 

where B, T], vp, I are the magnetic induction, viscosity, Fermi velocity and mean path length 
of the electron correspondingly, n m = B/<f) ~ 10 10 cm -2 one obtains for characteristics units 
of time, space and diffusion constant Kd(see EqJTOl) 



t — 5 • 10 8 sec, AL ~ 10 4 cm, n d = k ~ « (45) 

to 

the characteristic size of the interface in the dimension units L c ~ AL// 1 / 3 . 

For BCS superconductor where A — > the dimensional heating parameter Sa — 
n m v F /T 2 ~ i 2 n m < 1 (here C p ~ mp F T, while Bp and pp are the Fermi en- 
ergy and the momentum correspondingly. At low temperature T < Aq , where 
C p ~ (rappA^ 2 /T^ 2 J exp(— A /T) the heating parameter Sa grows dramatically Sa — 

^(^) 1/2 exp(A /T)»l. 

The heat parameter Sa is responsible for type of the instability. In particular at low tem- 
peratures (T << T c , where T c is the critical temperature) the parameter Sa is large and the 
instability develops on Mullens-Sekerka scenario (see Fig. 3) typical for dendritic instability 
{4]. On the other hand, at temperatures close to the critical, when the heating parameter 



Sa < 1 , the instability emerges as a periodic pattern (see figs. 4-6) 24]. 

Vortices and antivortices in the unstable pattern move with velocities (see Eq. (123]) ) 
u ~ 10 5 cm/ sec .If the difference of the "rapid" and "slow parts" of the front 5v ~ lOu 
the the flux pattern might reach the microscopic magnitudes L for very short time of the 
instability development r ~ 30/i sec, L w 3cm. The heat fluctuation for this time cannot 
significantly relaxes because it moves along the front on distance Sy ~ yfttr ~ 0.01cm. 

The main results of this paper are presented in Figs 3-6 where the increments of the insta- 
bility were drown for various anisotropy parameter a, heat conductivities inside and across 
the sample (k and r correspondingly) and heat annihilation coefficient Sa- We conclude that 
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heat released by the fTux-antifTux annihilation results in instability of the interface separating 
ffuxons and antifluxons areas even in the case of weak anisotropy of the superconductor. 

On the other hand if the superconductor is strongly anisotropic, the instability emerges 
even for weak heat. In more experimentally common case of moderated heat and anisotropy, 
both these mechanisms work together creating the instability of the flux-antiflux front. 

Our major conclusion is that the anisotropy of the superconducting layered structure 
alone cannot explain instability of the flux antiflux interface in weakly anisotropic materials 
as Nb and MgB 2 . From this point of view without the heating, the flux-antiflux front in the 
Nb superconductor (a ~ 0.9) without heating should be stable for any angle $ while strong 
heating destroys the front . If the heating caused by the vortex-antivortex annihilation is 
large then the vortex antivortex front instability should be detected in completely isotropic 
superconductors like MgB 2 (Fig. 6, curve 1). In fact it should be noted that even small 
heating might be essentially important to cause the instability. 

This theory is appropriate in the flux flow regime. The spatial disorder might affect the 
result by two different ways. It can both modify the linear profile of the magnetic induction 
at the front and affect the mechanism of the heat at the interface. 
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VII. APPENDIX I 



Substituting Eqs. ( 129]) -( 134)) in the initial set of Eqs. (I24l)- (|26j) . and use the stationary- 
solution in the form Eqs. ( lT4"|) - (fTo'j) one obtains for perturbations 



dt dx dx x dx dy 



v *djj =0 



dt dx 
86 



dSv x dip 



85v, 



'd 2 e d 2 9 s 



dt \dx 2 ' dy 2 



2S A (N (-b o iP)v o x + r9 = 



(46) 

(47) 
(48) 



Looking for the solution of these equations in the form Eq. (l32|) one obtains near the 
interface where b = —I 2 ^ 3 x, N = J 1 / 3 (see Eqs. fl20l) - fT211 ) the set of differential equations 
with uniform coefficients: 



A - J yy lh 2 ) 4> (x) + ( lxy + lyx ) ikl^^^- + lxJ la ^^- - lyjklk (x) (49) 

d( (x) 



— ry J 3 . 

/ XX 



86 (x) 

- ly jkie (x) - lxx i^M = o 



-2 lxy ikI^{x)-2 lxx I^ 
dip (x 



8x lyx 
dip (x 



dx 



dx 



+ (\ + 2 lxx I)((x)+ 7xx lle(x)=0 



(50) 



lxy ikI-*S A ^ ( x ) +lx jls A ^^-2 lxx IS A ( (x)+(A + r + Kk 2 - lxx I-*S A ) 6 (x)- K ^l = 



dx 2 
(51) 



These equations should be completed by the boundary conditions 



(52) 



x=—x c /2,x c /2 



The functions ip (x) , ( (x) and 9 (x) should be symmetrical and localized at the flux- 
antiflux interface where x < x c while x c ^> 1 is the cutoff where these functions go to zero 
(see Fig. 2 and Eq. flU). 

Looking for solution of Eqs. f l46|) -f l5T|) in the form 



c 

vV 



( (A n sinp n x + B n cosp n x) * 

(C n sinp n x + D n cosp n x) 
\ (E n sinp n x + F n cospnx) J 



(53) 
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one obtains equations for A n , B n ,C n} D n , E n , F n coefficients 



A 



C n 

Da 

E n 
\F n J 



(54) 



where matrix A reads 



— t\ 2 



p n 2j xy ikl 



-2j xy ikU 



2j xy i Pn kl3 -jyJkU Vnlxx 1 '" 



— 2 

^Vnlxx 1 ^ 



1xy ikIsS A 



Ixx^SaPu 



2 lxy ikli 



A + 2 lxx I 
A + 2 lxx I 



-Ixx^SPn 

^ xy ikllS A 



- 2 1xx IS A 







-lyx ikI 

-Ixx^n 

T ± 

'Y 1 3 

/ XX 



A + r + nk 2 — 
-1xx I±s Sa + kj% 



Vnlxx*- 

-lyx ikI 


I XX 



\ 



A + r + nk 2 — 



-2 lxx IS A 

-1xx^Sa + kpZ J 
(55) 

where p n = (2n + 1) vr/x c is obtained from the boundary condition ( |52l) . The dangerous 



harmonic with n = is responsible for instability. In this case po ~ x c 1 ~ ~ ^ f° r 
small slope of the magnetic induction at the interface (/< 1) this matrix can be simplified: 





\-^ y yIik 2 -J yx ikl3 

o x- lyy i\k 2 

-2 lxy ikll X + 2 lxx I 

-2 lxy ikll A + 2 7xx ./ 

2 

^ xy ikIiSA —2j xx IS A Atrtw-7 II i3,-i | 

l xy ikllS A -2 7aKB J5U 



-Jyx ikl3 



-lyx^I 


T J 3 



A + r + k/c 2 — 7 if 5/ 






ri 

T J 3 

/ XX 



A + r + nk 2 — ^ XX I^S/ 
(56) 



and can be represented as generated Jordan matrix. 
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K 
K 



where K is the 3x3 matrix 



K 



(A- 7to /^ 2 ) (-lyjkll^ {-l yx ikl) 
(-2 lxy ikl^j (X + 2 lxx I) (j x J%) 
(j xy ikI%S A ) (-2j xx IS A ) (\ + r + Kk 2 --f xx I^S / 



giving the equation for in the form 



a 3 + a 2 + n^ 2 ) + a (r 2 + u 2 k 2 + u 3 k 4 ) + (n 4 k 2 + n 5 fc 4 ) 



T 2 = 2 7xx ./r 

n 2 = 2 (-fl y ~ l XX lyy) it + 2 lx jK ~ lyyl^T + I^S A ^Jyy ~ ^ 

n 4 = 2 (-f 2 X y - J xx J yy ) rll + J xx J 2 xy S A I% 

n 5 = 2 (-fly - l XX lyy) 



VIII. APPENDIX II 



The solutions of the Eq. (|59|) at small k — > read: 

A = ~^k 2 + A k* 

1 2 

n 2 n 4 r^ 2 n 5 
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Ai 
Ci 

B l 



and 



1 



i\ + Jr?-4r 2 



ri - Vr? - 


■4r 2 




- n 2 ] + 2 [n 4 - r 2 ni] 






-4r 2 




+ 4r 2 





3Vr?-4r 2 -r! 


+ 


vn- 


-4r 2 




[25^! + n 3 ] + 2 [n 5 + Sxnx] 


4r 2 + ri 




-4r 





(61) 
(62) 

(63) 
(64) 



A 2 
C 2 = 

B 2 = 
Ao = 



C 2 + B 2 k 2 + A 2 k 4 ; 
1 



(65) 



-r 1 -,/r?-4r 2 



r\ + y/ri - 4r 2 


[n 2 - r^] + 2 [r 2 n x - n 4 ] 




y/rl - 4r 2 + r\ 


-4r 2 





- 4r 2 + Ti 


+ 


Vr? - 4r 2 + r\ 


[2B 2 n! + n 3 ] - 2 [n 5 + b 2 il 1 ] 




^V\-AT 2 + T 1 


-4r 2 



(66) 



(67) 
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Figure Captions 

Fig.l Geometry of the problem, i? is the angle between the x axis and the ab plain 
of the layered structure and c axis is perpendicular to the layers of the superconductor. 
Fig.2 Structure of the vortex-antivortex interface. 

Fig.3 Mullins-Sekerka instability for super large heat. The increment Re A versus k 2 (k 
is the wave vector along the ffux-antiflux front). 

Fig.4 The increment Re A versus k 2 for isotropic superconductors 
(a = 0.9, •§ = 7r/4)with different in-plane diffusion constant re. Curves 1,2,3 correspond 
to re — 0.1; 0.5; 1 respectively. Here 7 XX = 1 yy — 0.545, ^ xy = 0.055,/ = 0.5, r = 
0.0148, Sa = 0.89. Instability disappears as the diffusion constant grows. 

Fig. 5 The increment Re A versus k 2 for anisotropic superconductor with different re- 
laxation constant r . Curves 1,2,3 correspond to parameters r = 0.049; 0.13; 2.66 respectively. 
The instability disappears as the relaxation parameter grows. Here j xx = 7 y?/ = 0.545, ^ xy = 
0.055, re = 0.1, 1 = 0.5, S A = 0.89. 

Fig.6 The increment of instability Re A versus k 2 for different anisotropy. Curve 1 
for isotropic superconductor (a = 0.9, $ = 7r/4, ^ xx = ^ yy = 0.545, ^ xy = 0.055) shows 
the instability at heating coefficient Sa = 0.89 while curve 2 demonstrates the lack of 
instability at heating coefficient Sa = 0.8. The curve 3 exhibits instability for anisotropic 
superconductor (a = 0.1, i? = 7r/4, ^ xx = = 0.55, -y xy = 0.45) and even more small 
heating coefficient S A = 0.1. (Here r = 0.0148, re = 0.1, / = 0.5). 

Fig. 7 Qualitative picture of the Kelvin- Helmholtz instability at the vortex-antivortex 
interface. 
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This figure "Figl.jpg" is available in "jpg" format from: 
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